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ROBUSTIFYING  THE  KALMAN  FILTER; 


Protection  Against  Symmetrically  Straggling  Measurement  Errors 

D.  P.  Gaver 
P.  A.  Jacobs 


1.  INTRODUCTION. 

Tracking  and  prediction  algorithms  based  on  simple  Gaussian  (normal 
distribution)  measurement  errors  and  structural  models  are  commonly  used  in 
practice  under  the  name  of  KALMAN  Filters.  If  (a)  measurement  errors  are  not 
suitably  Gaussian,  e.g.,  if  occasional  outliers  occur  or  (b)  true  structural 
behavior  is  not  simple,  perhaps  displaying  apparently  discontinuous  behavior 
caused  by  unfavorable  sensor-target  orientation,  then  traditional  Filter 
performance  may  dramatically  degrade.  In  this  paper,  we  will  propose  and 
study  procedures  based  on  an  elaborated  model  of  the  KALMAN-type  but  with 
the  measurement  errors  coming  from  a  family  of  possibly  suitable  non-Gaussian 
distributions  (e.g.,  Student-t)  to  represent,  and  suitably  compensate  for  more- 
thick-tailed-than-Gaussian  measurement  error,  i.e.,  distributions  with  long 
straggling  tails  having  the  tendency  to  produce  symmetric  outliers. 

In  particular  the  basic  stochastic  model  considered  here  is 


—  ®n-l  +  (LI) 

Yn  =  ®n  +  En  0  -2) 

where  {(On)  are  independent  normal/Gaussian  random  variables  with  mean  0 
and  variances  { xn )  and  {£„}  are  independent  random  variables  having  mean  0. 
The  random  variable  9n  is  unobservable.  The  random  variable  Yn  is  interpreted 
as  the  observation  of  0n  made  with  measurement  error  en;  £n  is  not  Gaussian, 
but  controllably  long-tailed.  The  problem  is  to  estimate  9n  from  Yi,...,  Yn  in 
the  simple  recursive  fashion  that  characterizes  the  classical  KALMAN  filter. 
Expression  (1.1)  is  a  simple  random  walk  and  does  not  represent  very 
interesting  dynamics,  but  does  provide  suggestive  illustrations. 


In  the  next,  or  second,  section,  we  will  describe  a  procedure,  the  ALMA 
(standing  for  KALMAN  with  outliers  suppressed),  which  is  based  on  a  model  in 
which  the  components  of  the  error  sequence  (en)  have  a  Student-t  distribution. 


7»vTaT'C^xr^T« TXrv'v 


In  the  third  section,  the  traditional  KALMAN  procedure  will  be 
described.  It  is  based  on  the  assumption  that  components  of  (en)  have  iid 
normal  distributions.  Finally,  a  robust  procedure  due  to  West  [  1981  ]  will  be 
described. 


In  section  4  results  of  an  extensive  simulation  experiment  will  be 
presented  and  discussed.  The  simulation  experiment  compares  the  various 
procedures.  The  results  indicate  that  the  ALMA  procedure  is  significantly 
better  than  the  KALMAN  when  the  true  measurement  error  distribution  is 
Student-t.  Further,  there  is  not  much  lost  in  using  the  ALMA  procedure  instead 
of  the  KALMAN  when  the  true  measurement  error  distribution  is  normal. 

2.  THE  ALMA  FILTER  AND  RELATED  PROCEDURES. 


While  many  measurement  errors  of  physical  quantities  are  approximately 
normal,  especially  "in  the  middle"  of  their  distribution,  there  can  well  be 
thicker-than-normal/Gauss  tails  and  also  occasional  extreme  outliers;  that  these 
can  have  seriously  degrading  effects  in  regression-like  problems  has  been  the 
subject  of  considerable  research;  we  cite  books  by  Mosteller  and  Tukey  (1977), 
Huber  (1981),  Hampel  (1986);  in  the  time-series  context  the  article  by  Martin  and 
Yohai  (1986),  which  contains  many  references;  also  lately  the  articles  by  West 
and  his  associates  (1981,1985);  it  is  to  West's  approach  that  our  methodology 
should  best  be  compared. 


One  way  to  model  these  features  is  to  extend  the  tails  of  the  normal  by 
continuous  scale  mixing.  Such  an  approach  can  lead  to  the  Student-t  form,  and 

to  many  other  useful  forms  as  well.  We  will  assume  here  that  { en }  are 
independent  random  variables,  now  having  in  the  Student-t  distribution  with 
mean  0,  scale  Gn  {not  the  standard  deviation)  and  d  degrees  of  freedom;  that  is. 


pe  (u)  =  c(d)  — 
n  on 


(2.1) 


Let  y,  denote  the  i111  measurement  and  yn  =  (yj . yn).  Assume  that  0n- 

ll  yn-*  has  a  normal  distribution  with  mean  mn_i  and  variance  Cn-i.  Since  con  is 
assumed  to  have  a  normal  distribution  with  variance  Tn,  0nl  yn' 1  has  a  normal 


distribution  with  mean  mn-i  and  variance  Cn=Cn-i+Tn.  Thus,  from  (1.1),  (1.2), 
and  (2.1) 

P  {0ne d0,Yne dy  I  Y1=y1,...,Yn-i=yn.i ) 


■xp |t '8“n#l)  -i(d+1)ln  1+(^f)  ?  d0dy 


(2.2) 


“  K  exP  [4 (9  ^(y)^  +  jQ(y)]d0dy 


where  the  approximation  replaces  the  expression  in  the  exponent  by  an 
approximating  quadratic  in  0. 


The  ALMA  procedure  provides  a  Gaussian  approximation  to  the 
distribution  of  0n  |  yn,  but  one  that  emphatically  differs  from  the  classical  linear- 
in-observations  form.  Following  an  argument  in  Gaver  et  al.  [1986], 
differentiate  both  sides  of  (2.2)  with  respect  to  0  to  obtain 


e-p(y)  0-mn. !  d+1  0-y 

C(y)  n#  +  ~3~  2 


0-yf  l 


(2.3) 


Equating  the  terms  involving  0  results  in  the  following  equation: 
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where  the  weight 


Furthermore,  equating  the  constant  terms  results  m 


The  ALMA  procedure  approximates  0n  I  yn  by  the  normal  distribution 
having  mean 


K-ipr  +  w(yn,yn4 
n  n 


and  variance 


Cn  =  [-“r  +  w(yn)  4r 


where 


w(yn)  = 


d-H  1 


Note  that  the  weight  w(yn)  involves  the  unknown  0.  One  implementation 
uses  approximate  weights  of  the  form 


wk(yn)  - 


l+p^ 


™n-Q2  k 


(2.1C 


When  k=l,  mn.i  is  used  in  place  of  0  in  (2.9). 
When  k=j,  0.5(mn_i  +yn)  is  used  in  place  of  0. 
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The  basic  ALMA  procedure  is  to  evaluate  w^yn)  and  then  use  it  to  find 


1  .  ,  1  -) 
73  *  Wk<yn>  -J 

'"'n  °n 


(2.11) 

(2.12) 


The  point  estimate  of  0n  given  yn  is  0n=mn  and  m  estimate  of  the  variance 

of  0n  is  Cn.  Thus  the  procedure  provides  a  particular  Gaussian  posterior 
approximation.  In  other  similar  contexts,  non-linear  filters  for  example,  it  has 
been  suggested  that  the  procedure  (2.10)  -  (2.12)  be  iterated  with  the  newly- 
computed  mn,  replacing  mn_i  in  (2.10)  -  (2.12)  in  each  iteration.  In  the 
simulations  0, 1  and  2  iterations  were  implemented,  and  the  results  compared. 


The  ALMA  procedure  is  an  iterative  reweighting  procedure.  In  the 
ordinary  regression  context  another  weight  has  been  suggested:  the  so-called 
(Tukey)  biweight,  cf.  Mosteller  and  Tukey  (1977).  In  our  context,  the  biweight 
procedure  can  replace  the  weight  w^(y)  in  the  ALMA  procedure  with  the 
biweight 


(2.13) 


The  variance  of  a  Student-t  distribution  with  d  degrees  of  freedom  and 
scale  a  is  d>3,  otherwise  being  infinite.  Hence  the  (bi)weight  wB(y)  uses 

the  measurement  y  if  I  y  I  is  within  a  standard  deviations  of  mn_i,  the  estimate  of 
0n.l.  The  weight  is  zero  if  the  deviation  is  greater. 

As  was  done  in  the  basic  ALMA  procedure,  0,  1,  and  2  iterations  of 
(2.10)-(2.12)  were  tried,  with  wB(yn)  replacing  wic(yn),  for  values  of  a=5,  7,9 
and  k=l,  0.25. 


«.*  I.‘  . 


2.3  Aspects  of  the  Likelihood  Procedure. 

It  is  possible  for  the  likelihood  function  (2.2)  to  exhibit  two  local  0- 
maxima.  In  such  a  case,  the  likelihood  procedure  approximates  the  local 
maxima  and  chooses  the  one  which  globally  maximizes  the  likelihood. 

To  examine  the  details  let 

f(0)  =  -^-lnP{e  ed6,Y  edy} 

HA  1  n  n  31 


f 

-(e-mn-i) 

< 


d  + 1 

+  T" 


1  + 

L  V 


<*n 


-1 


A 


y_± 


d9dy. 


J 


(2.14) 


Now  it  is  clearly  possible  for  f(0)=O  to  have  multiple  roots.  To  be 
specific,  f(0)=O  for  those  0  satisfying 


O  =  03+e2(-2y-rn  ) 


(2.15) 


+  - 


®[^d  +  y2+W+l)C;  +  2yaVlj 

[•mn.i°5d-ni,.iy2-<d+,»'c;} 


The  properties  of  this  cubic-in-0  equation  can  be  deduced  from  classical 
results. 


i ^  v, a  .  «*,  ■»*„  ^ »*,  «*. v,  •’ . 


2 


Cn  21 

2d2-5d(d+l)  — r  --Hd+l) 
erf 


* 


-i  3 


d+(d+l) 

°n 


then  if 


D>0  (2.15)  has  1  real  root  and  two  conjugate  imaginary  roots; 
D=0  (2.15)  has  3  real  roots,  at  least  two  of  which  are  equal; 
D<0  (2.15)  possesses  3  real  and  unequal  roots. 


Note  that  if  d=°°  so  that  en  has  a  normal  distribution,  then  certainly  D>0 
and  (2.2)  has  a  unique  maximum.  If  d«»  and  C*on2  is  small  enough,  then  D>0 

and  once  again  (2.2)  will  have  a  unique  maximum.  If  d<«>  and  C^an2  is  large 


enough  (actually,  larger  than 
2 

values  of  (y-mn.,)  and  (2.15) 
have  two  local  maxima. 


(6V3-10)d  4 d  x 

d+1  ~  d+r’ 


then  D<0  for  an  interval  of 


will  have  3  real  unequal  roots;  in  this  case  (2.2)  will 


The  likelihood  procedure  computes  D.  If  D>0  it  uses  the  ALMA 
procedure  with  weight 


-l 


,  ,  d+1 
wk(y)»-g- 


(2.17) 


to  compute  0n.  If  D<0,  then  two  candidate  estimates  0b  and  02  of  0  are 
computed.  Both  estimates  are  obtained  via  the  ALMA  procedure  (2.7)-(2.9). 
One  approximates  weight  (2.9)  by  setting  0=010.!  as  in  (2.10);  think  of  the  result 
as  prior-dominated.  The  other  approximates  weight  (2.9)  by  setting  0=y,so  that 

w(y)  =  ;  the  result  is  data-determinated.  The  likelihood  function  is  then 

evaluated  at  each  value  of  0:0!  and  02  The  quoted  estimate  of  0n  is  set  equal  to 
the  0j  that  comes  closest  to  maximizing  the  global  likelihood;  the  estimate  of  the 
variance  is  set  equal  to  the  corresponding  Cn. 

3.  THE  KALMAN  AND  WEST  PROCEDURES. 

In  this  subsection,  the  traditional  KALMAN  procedure  will  be  described 
for  the  model  (1.1)-(1.2).  A  procedure  proposed  by  West  (1981)  will  also  be 
discussed. 

3.1  The  KALMAN  Procedure. 

The  KALMAN  filter  finds  the  estimate  6n  of  0n  which  minimizes  the 
conditional  mean  square  error  of  (6n'9n)  given  yn.  If  { en }  are  independently 
normally  distributed  with  mean  0  and  variances  {yn}»  then  the  KALMAN  filter 
can  be  viewed  as  a  Bayesian  updating  procedure;  see  Meinhold  and  Singpurwalla 
(1983). 

The  Bayesian  KALMAN  procedure  assumes  0n.!lyn  l  is  normal  with  mean 
nrin.i  and  variance  C,,.!.  Thus,  from  (1.1)  0nlyn  l  is  normal  with  mean  mn_,  and 
variance  C*=Cn_i  +xn.  From  (1.2) 


2  2" 

P(0ne  d0,YnG  dy  I  yn  l)  =  K  exp  \  ^  — —  -  d0dy 


(3.1) 


West  proposes  an  estimation  procedure  for  0n  given  yn  in  the  case  in 
which  the  density  p^  is  symmetric  about  0.  In  the  special  case  in  which  p^  is 

normal.  West's  procedure  reduces  to  the  KALMAN  filter. 


Once  again,  assume  0n_i  yn_1  is  normal  with  mean  mn_i  and  variance 


Cn_i  so  that  ©nly0'1  is  normal  with  mean  inland  variance  Cn=Cn_i+  xn 


P{0  ed0,  Ynedy|Yi=yi . Yn.i=yn.i) 


=  Kcxp  -r(0-mn.if  — r  +  ln  p  (y-0) 
z  r  n 

L  '■'n 


(3.6) 


(9mn-i)" 


K  exp  ^0-mn.i)  —  +1  In  P^Cy-mn-l)  +  g(y-mn_i)(0-mn.i)  -  G(y-mn„i) - - - I  d0dy  (3.7) 


where  a  Taylor  expansion  provides 

.  -d 


g(u)  =  du  \ 


(3.8) 


G(u)  =  ^%n(u)‘ 


(3.9) 


Completing  the  square  in  (3.7)  results  in 


P  [0ne  d0,Yne  dy  I  Y 1  =y  1 Yn- 1  =yn- 1 1 


ir  i  v  r  1  T  V 

«Rexp^  —  +G(y-mn-i)  (9-mn-i)-g(y-mn.i)  — +  G(y-mn-i)  L  (3.10) 

Cn  Gn 


Hence,  P{0n€d0  I  Y,=y|, 
distribution  having  mean 


.Yn=yn )  is  approximated  by  a  normal 


mn  —  mn-i  +  Cng(yn-mn.| ) 


(3.11) 


*  0  V  ^  C  C  «%/  «/V  O  4 *  *  *  V  O-T.* 


and  variance 


-l 


Cn  = 


+ 


G(yn-mn-l) 


(3.12) 


In  the  special  case  in  which  en  has  a  Student-t  distribution  with  d  degrees 
of  freedom  and  scale  parameter  on. 


(u)=c(d)-~- 

On 


-(d+1) 

ll  2 


(3.13) 


and 


2  ' 

-2  - 

»  d+1  1 

0(u)  =  -j-a2 

h?)* 

(3.14) 


(3.15) 


Since  G(yn-mn-i)  is  playing  the  role  of  a  variance  in  (3.10),  but  may 
become  embarrassingly  negative  for  large  u.  West  suggests  that  it  be  replaced  by 
max(0,G(yn-mn„i));  this  step  has  been  taken  in  the  simulations  that  illustrate  the 
various  procedures  proposed  here.  West  suggests  another  possibility  in  West  et 
al.  (1985). 

4.  A  SIMULATION  EXPERIMENT. 

All  simulations  were  carried  out  on  an  IBM  3033AP  computer  at  the 
Naval  Postgraduate  School.  Random  numbers  were  generated  using  the 
LLRANDOMII  random  number  package;  cf.  Lewis  and  Uribe  (1981). 

For  each  replication  of  the  simulation  the  model  of  (1.1)-(1.2)  is  generated 
for  n=0,l,...,100.  In  the  simulations  reported  below  { con }  are  iid  normal  with 


mean  zero  and  variance  one.  For  each  replication,  estimates  &n  of  0n  given  yn 
are  computed  using  each  of  the  procedures  described  above.  The  data  collected 
are  the  estimation  error  &n  -  0n  for  n=25,  50, 75,  100  and  the  estimate  of 
variance  Cn,  n=25,  50,  75,  100.  The  number  of  independent  replications  is 
1000. 

Tables  1  and  2  report  results  of  the  KALMAN  and  ALMA  procedures  for 
simulations  in  which  {en}  are  iid  normal  with  mean  zero  and  variance  one.  The 
ALMA  procedure  actually  uses  the  incorrect  measurement  error  model  that 
{en}  are  iid  Student-t  with  d=3  degrees  of  freedom  and  variance  equal  to  one. 
Results  for  the  ALMA  procedure  are  shown  for  weights  as  in  (2.10),  for  k=1.0 
and  k=0.25.  The  procedure  was  iterated  0,  1 ,  and  2  times. 

Table  1  shows  statistics  of  6n-0n  for  n=25,  50,  75,  100.  As  anticipated,  the 
KALMAN  procedure  which  uses  the  correct  (normal)  model  exhibits  the 
smallest  variance  of  6n  -  0n.  The  ALMA  procedure  with  k=0.25  and  0  iterations 
and  the  ALMA  procedure  with  k=l  and  1  iteration  have  the  smallest  variances 
for  the  ALMA  procedures. 


Table  2  exhibits  the  estimates  of  the  variance  of  0n,  namely  Cn,  for  the 
ALMA  procedure  for  n=25,  50,  75,  100.  The  KALMAN  estimate  of  the 
variance  is  the  constant  0.618  for  all  of  these  n.  This  constant  is  the  limiting 

lim 

solution  to  equation  (3.4)  with  Tn  =  Yn=l  >  ^at  *s»  C=  Cn 


n-  >°° 


a  simple  quadratic  with  appropriate  solution 


1+V5 

C  =  o.618. 


The  variance  of  6n-0n  for  the  KALMAN  procedure  in  Table  1  is  close  to 
the  calculated  0.618. 


The  mean  values  of  Cn  for  the  ALMA  procedure  with  k=0.25  and  0 
iterations  and  k=l  with  1  iteration  are  about  half  that  of  the  corresponding 
variances  of&n-0n  in  Table  1. 


Tables  3-4  report  results  for  a  simulation  in  which  (en)  are  iid  Student-t 
with  3  degrees  of  freedom  and  variance  equal  to  1.  Table  3  reports  statistics  of 
the  estimation  error.  6n  0n,  for  t  ie  KALMAN,  ALMA,  Biweight,  Likelihood, 
and  West  procedures.  As  usual,  the  KALMAN  procedure  assumes  { en }  are  iid 
normal  w  ith  mean  0  and  variance  1 .  The  other  procedures  assume  { en }  are  iid 
Student-t  with  3  degrees  of  freedom  and  variance  equal  to  1 .  The  ALMA 
procedure  with  k-  33  and  no  iterations  exhibits  the  smallest  variance  of  §n  -  0n. 
The  more  complicated  Likedhood  procedure  with  k-0.25  and  no  iterations 
exhibits  the  next-smallest  variance.  The  ALMA  with  k=l  and  1  iteration 
exhibits  the  third  smaUeo  . unan.;e. 

The  Biweig.it  procedure  was  implemented  with  the  constants  in  the  weight 
(2.13)  a=5,7.9  and  k~.23  and  1.  the  procedure  was  iterated  0,1,  and  2  times. 

The  results  for  a=-5  were  much  worse  than  those  for  a=7  and  9  indicating  that 
a=5  is  not  large  enough  to  suppress  outlying  values;  they  are  not  reported. 
Iterating  the  biweight  procedure  I  and  2  times  did  not  improve  the  results  for 
any  values  of  a.  'Hie  results  of 'fable  7  indicate  that  the  biweight  procedure  with 
the  smallest  variance  uses  k--l  .0  and  a=7  with  no  iterations. 

7 Tie  West  procedure  deescribed  in  West  (1981)  as  currently  implemented 
does  not  do  as  well  as  the  KALMAN.  The  statistics  of  Cn  in  'fable  4  seem  to 
indicate  that  the  difficulty  is  with  the  estimate  of  variance,  Cn;  the  fix  for 
negative  G(y-mn  , )  makes  it  possible  for  Cn  to  increase  by  one  in  successive 
times  over  long  periods  of  time. 

Table  4  exhibits  the  statistics  of  Cn.  The  KALMAN  procedure,  the 
ALMA  procedure  with  k=0.25  and  (J  iterations,  the  ALMA  procedure  with  k=l 
and  1  iteration,  the  Likelihood  procedure  with  k=0.25  and  0  iterations  and  the 
Biweight  with  k-  1 ,  a=7  all  have  mean  Cn  approximately  half  the  variance  of  §n- 

0n. 

5.  CONCLUSIONS. 

'fhe  simulation  reviles  >  burned  to  date  indicate  that  a  satisfactory  robust 
KALMAN-AI.MA  picretiiiii'  utilizes  the  k-0.25  weight-starting  option  and 
requires  no  iteration.  Wiui<-  due  above  tiltei  in  about  7',’f  less  eft icient  than  the 
**  KALMAN  when  measurement  errors  are  ideally  Gaussian,  it  is  about  6%  more 

efficient  when  enors  are  long  aided  non  Gaussian;  efficiency  is  in  terms  of 
ratios  of  ''estimated  )  variances  and  :s  not  the  only  meaningful  criterion. 


.-■v v.v.; 


Examination  of  Table  3  reveals  through  values  of  skewness,  and  kurtosis  ,  that  as 
anticipated,  the  robust  ALMA  estimation  errors  are  substantially  more  closely 
Gaussian  than  are  the  corresponding  KALMAN  products  when  measurement 
errors  are  Student-t. 
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Table  1 

Statistics  of  §n-0n 

Normal  Measurement  Errors  with  Variance  1 


Time  n: 


Proc  Nbr  k 
Iter 


25 


M  V  S  K 

0.00  0.61  0.04  0.01 
0.02  0.91  0.07  0.60 
0.00  0.65  0.02  0.07 
0.01  0.70  0.05  0.02 
0.02  0.70  0.00  0.07 
0.01  0.71  0.02  0.09 
0.02  0.77  0.01  0.11 


50  75  100 


MVS  KMVS  KMVS  K 


0.02  0.62 
0.03  0.77 
0.03  0.65 
0.04  0.69 
0.03  0.76 
0.04  0.75 
0.02  0.83 


0.01  0.15 
0.21  1.40 
0.04  0.03 
0.15  0.62 
0.02  0.09 
0.06  0.06 
0.02  0.10 


0.03  0.63 
0.01  0.78 
0.02  0.64 
0.01  0.68 
0.03  0.74 
0.02  0.72 
0.03  0.82 


0.02  0.29 
0.01  0.18 
0.01  0.14 
0.00  0.02 
0.04  0.10 
0.01  0.04 
0.04  0.09 


0.00  0.60 
0.01  0.84 
0.00  0.67 
0.03  0.73 
0.00  0.75 
0.00  0.76 
0.00  0.81 


-0.05  -0.19 
-0.19  0.78 
-0.05  41.17 
-0.06  -0.01 
0.02-030 
-0.01  -0.25 
0.01-031 


K  =  KALMAN 
A  =  ALMA 


Mean 

Variance 

Skewness 

Kurtosis 
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Table  2 
Statistics  of  Cn 

Normal  Measurement  Errors  with  Variance  1 


Time  n 

25 

50 

75 

100 

Proc 

Nbr 

k 

M 

V 

M 

V 

M 

V 

M 

V 

Iter 

A 

0 

1.0 

.50 

.08 

.49 

.07 

48 

.06 

.48 

07 

0.25 

.31 

.01 

.30 

.01 

.30 

.01 

.30 

01 

l 

1.0 

.23 

.02 

.22 

.02 

.22 

.02 

.22 

02 

0.25 

.04 

.00 

.14 

.00 

.14 

.00 

.13 

00 

2 

1.0 

.14 

.01 

.13 

.01 

.13 

.01 

.13 

01 

0.25 

.09 

.00 

.09 

.00 

.09 

.00 

.09 

00 

Procedure  (Proc.)  Statistics 

A = ALMA  M=Mean 

V  =  Variance 
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Table  3 

Statistics  of  §n-0n 

Student-t  Measurement  Errors  with  3  degrees  of  freedom  and  Variance  1. 


Time  n: 

Proc  Nbr 
Iter 
K 

A  0 


_ 25 _ 

50 

75 

100 

k  a 

M 

V 

S  K 

M  V  S  K 

M  V  S  K 

M  V  S  K 

. 

0.01 

0.57 

0.48  17 

0.02  053  0.02  22 

0.02  0.67  0.78  8.7 

0.04  0.54  0.17  1.7 

1.0  - 

0.03 

0.67 

0.07  1.0 

0.02  058  0.08  1.1 

0.02  0.71  0.02  1.4 

0.01  0.65  0.10  1.7 

0.25  - 

0.01 

0.53 

0.16  1.6 

0.01  0.48  0.08  15 

0.02  057  0.01  21 

0.02  0.50  0.26  23 

1.0 

0.01 

0.55 

0.09  1.5 

0.01  0.49  0.09  1.2 

0.02  0.61  0.01  1.8 

0.01  0.52  0.23  24 

0.25  • 

-0.01 

0.63 

0.46  4.1 

0.01  058  0.03  3.1 

0.03  0.69  0.06  35 

0.03  0.58  0.20  2.7 

1.0 

-0.01 

0.58 

0.16  2.4 

0.01  054  0.05  1.6 

0.03  0.64  0.02  1.9 

0.02  0.54  0.26  2.6 

0.25  - 

-0.02 

0.71 

0.66  5.5 

0.01  066  0.09  43 

003  0  79  0.18  4  .9 

0.04  0.64  0.05  2.9 

1.0  7 

0.01 

0.56 

0.17  15 

0.01  0.51  0.07  22 

002  059  0.09  2.6 

0.02  0.52  0.31  24 

0.25  7 

0.01 

0.61 

0.68  5.2 

0.01  057  0.04  3.9 

0.03  0.69  035  6.0 

004  0.56  0.13  2.7 

1.0  9 

-0.01 

0.57 

0.42  3.5 

0.01  054  0.00  3.0 

0.02  0.62  0.07  33 

0.03  0.53  0.22  25 

0.25  9 

-0.01 

0.63 

0.76  5.7 

0.01  057  0.05  4.1 

0.03  0.72  0.75  9.9 

0.04  0.57  0.11  28 

1.0  ■ 

0.03 

0.68 

0.12  15 

0.03  0.56  0.09  1.0 

0.03  0.69  0.02  15 

0.01  0.65  0.13  1.8 

0.25  - 

0.01 

0-54 

0.23  1.7 

0.01  0.48  0.11  15 

0.03  0.60  0.04  21 

0.02  0.52  0.29  25 

1.0  • 

0.01 

0.55 

0.13  1.7 

0.01  0.49  0.09  1  2 

0.02  0.61  0.01  1.8 

0.01  0.52  0.23  24 

025 

-0.01 

063 

053  4  1 

0.00  057  0.07  28 

0.03  0.69  0.07  3.4 

0.03  0.58  0.20  27 

1.0 

0.00 

0.59 

022  25 

0.01  0.53  0.08  16 

003  053  0.02  1.9 

0.02  0.54  0.26  26 

025 

0.02 

0.72 

0.70  55 

0.00  0.65  0.01  3.9 

0.03  0.79  0.18  4.9 

0.04  0.64  0.05  29 

0.06 

112 

0.26  3.2 

051  384  0.14  43 

0.08  774  0.03  5.1 

0.79  1249  0.25  5.2 

K  =  KALMAN 
A  =  ALMA 
B  =  Biweight 
L  =  Likelihood 
W  =  West 


SiatiMics 
M  =  Mean 
V  =  Variance 
S  =  Skewness 
K  =  Kurtosis 


Table  4 
Statistics  of  Cn 

Student-t  Measurement  Errors  with  3  degrees  of  freedom  and  Variance  1 


Time  n: 

25 

50 

75 

100 

Proc 

Nbr  k 

a 

M 

V 

M 

V 

M 

V 

M 

V 

Iter 

A 

0  1.0 

- 

.46 

.07 

.44 

.06 

.48 

.07 

.45 

.06 

0.25 

- 

.29 

.01 

.28 

.01 

.29 

.01 

.29 

.01 

A 

1  1.0 

- 

.21 

.02 

.20 

.02 

.21 

.02 

.21 

.02 

0.25 

- 

.13 

.00 

.13 

.00 

.14 

.00 

.14 

.00 

A 

2  1.0 

- 

.13 

.01 

.12 

.01 

.13 

.01 

.13 

.01 

0.25 

- 

.09 

.00 

.09 

.00 

.09 

.00 

.09 

.00 

B 

0  1.0 

7 

.29 

.01 

.29 

.00 

.29 

.01 

.29 

.00 

0.25 

7 

.27 

.00 

.27 

.00 

.27 

.00 

.27 

.00 

B 

0  1.0 

9 

.28 

.00 

.28 

.00 

.28 

.00 

.28 

.00 

0.25 

9 

.27 

.00 

.27 

.00 

.27 

.00 

.27 

00 

L 

0  1.0 

- 

.44 

.06 

.44 

.06 

.46 

.06 

.46 

.06 

0.25 

- 

.29 

.01 

.28 

.01 

.29 

.01 

.29 

.01 

L 

1  1.0 

- 

.21 

.02 

.20 

.02 

.21 

.02 

.21 

.02 

0.25 

- 

.14 

.00 

.14 

.00 

.14 

.00 

.14 

.00 

L 

2  1.0 

- 

.13 

.01 

.12 

.01 

.13 

.01 

.13 

.01 

0.25 

- 

.09 

.00 

.09 

.00 

.09 

.00 

.09 

.00 

W 

-  - 

- 

8.8 

64 

16 

240 

23 

543 

29 

946 

Procedures  (Proc.) 
A = ALMA 
B=Biweight 
L= Likelihood 
W  =  West 


Safeties 
M=Mean 
V= Variance 
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